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Abstract 



In this paper, we consider a large class of purely discontinuous rotationally symmetric Levy 
processes. We establish sharp two-sided estimates for the transition densities of such processes 
^sO killed upon leaving an open set D. When D is a K-fat open set, the sharp two-sided estimates are 

given in terms of surviving probabilities and the global transition density of the Levy process, 
i— —i When D is a C 1,1 open set and the Levy exponent of the process is given by >J>(£) = 0(|£| 2 ) 

with cf> being a complete Bernstein function satisfying a mild growth condition at infinity, our 
two-sided estimates are explicit in terms of 'J, the distance function to the boundary of D and 
the jumping kernel of X, which give an affirmative answer to the conjecture posted in |16j . Our 
results are the first sharp two-sided Dirichlet heat kernel estimates for a large class of symmetric 
Levy processes with general Levy exponents. We also derive an explicit lower bound estimate 
'— 1 for symmetric Levy processes on M. d in terms of their Levy exponents. 
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5_i 1 Introduction 

a 

Due to their importance in theory and applications, fine potential theoretical properties of Levy 
processes have been under intense study recently. The transition density p(t, x, y) of a Levy process 
is the heat kernel of the generator of the process. However, the transition density (if it exists) of 
a general Levy process rarely admits an explicit expression. Thus obtaining sharp estimates on 
p(t, x, y) is a fundamental problem both in probability theory and in analysis. 
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The generator of a discontinuous Levy process is an integro-differential operator and so it is 
a non-local operator. Recently, quite a few people in PDE are interested in problems related to 
non-local operators; see, for example, (H [7] [M] and the references therein. 

When X is a symmetric diffusion on M rf whose infinitesimal generator is a uniformly elliptic and 
bounded divergence form operator, it is well-known that p(t, x, y) enjoys the celebrated Aronson's 
Gaussian type estimates. When X is a pure jump symmetric process on sharp estimates on 
p(t, x, y) have been studied in [81 [9], [10J [TTJ, [18] recently, which can be viewed as the counterpart of 
Aronson's estimates for non-local operators. 

Due to the complication near the boundary, two-sided estimates for the transition densities of 
discontinuous Levy processes killed upon leaving an open set D (equivalently, the Dirichlet heat 
kernels) have been established very recently for a few particular processes only. The first of such 
estimates is obtained in [11], where we succeeded in establishing sharp two-sided estimates for the 
heat kernel of the fractional Laplacian A a//2 := — (— A) a//2 with zero exterior condition on D c (or 
equivalently, the transition density of the killed symmetric a-stable process) in any C 1 ' 1 open set 
D. The approach developed in [11] provides a road map for establishing sharp two-sided heat kernel 
estimates of other jump processes in open subsets of R d . The ideas of [HJ were adapted to establish 
sharp two-sided heat kernel estimates of relativistic stable processes and mixed stable processes in 
C 1 ' 1 open subsets of M. d in |144 [T3] respectively. In all these cases, the characteristic exponents 
of these Levy processes admit explicit expressions, the boundary decay rates of the Dirichlet heat 
kernels are suitable powers of the distance to the boundary. On the other hand, a Varopoulos 
type two-sided Dirichlet heat kernel estimate of symmetric stable processes in «>fat open sets was 
derived in [1] ; this type of estimates is expressed in terms of surviving probabilities and the global 
transition density of the symmetric stable process. 

The objective of this paper is to establish sharp two-sided estimates on the transition density 
PD(t,x,y) for a large class of purely discontinuous rotationally symmetric Levy processes. Unlike 
the cases considered in [3] [HJ [13], [H] , the characteristic exponents of the symmetric Levy processes 
considered in this paper are quite general, satisfying only certain mild growth condition at oo. 
Moreover, the boundary decay rate of pn(t,x,y) is no longer some power of the distance to the 
boundary. The analysis of the precise boundary behavior of po(t,x,y) is quite challenging and 
delicate. The main tools to obtain the precise boundary behavior of pu(t,x,y) are two versions 
of the boundary Harnack principle obtained in [19} 121]. In this paper we combine the approaches 
developed in [U [UJ with these boundary Harnack principles to obtain sharp two-sided estimates 
for po(t, x, y), which cover the main results in [3] [TT] [TBI E] an d much more. 

Suppose that S = (St : t > 0) is a subordinator with Laplace exponent <p, that is, S is a 
nonnegative Levy process with Sq = and E [e _A5 ''] = e~*^M for every t, A > 0. The function (j> 
can be written in the form 



where b > and [i is a measure on (0,oo) satisfying A t)p,(dt) < oo. The constant b is called 

the drift of the subordinator and [i the Levy measure of the subordinator (or of <p). The function 
is a Bernstein function, i.e., it is C°°, positive and (— l) n-1 D n i^> > for all n > 1. In particular, 
since (f)(0) = and cp" < 0, the Bernstein function <p has the property that 






for all A > 1 and r > 0. 
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The Laplace exponent <fi is said to be a complete Bernstein function if the Levy measure /i of cj) 
has a completely monotone density fj,(t), i.e., (— > for every non- negative integer n. For 
basic results on complete Bernstein functions, we refer the reader to [23]. 

Throughout this paper, we assume that <p is a complete Bernstein function satisfying the fol- 
lowing growth condition at infinity (see |27|): 

(A): There exist constants 61,62 G (0, 1), a\ £ (0, 1), ai £ (l,oo) and Ro > such that 

aiA 5l </>(r) < ^(Ar) < a 2 A <52 ^(r) for A > 1 and r > i? . 

Note that it follows from the upper bound condition in (A) that <j) has no drift. 

Let W = (Wt : t > 0) be a Brownian motion in M. d independent of the subordinator S. The 
subordinate Brownian motion Y = (Yt : t > 0) is defined by Yf := Ws t , which is a rotationally 
symmetric Levy process with Levy exponent c/>(|£| 2 ). The infinitesimal generator of Y is C Y := 
—</>(— A). Here and below for a function ip on [0, 00), A) is defined as a pseudo differential 
operator in terms of Fourier transform; that is, ip(—A)f(£) := —ip(\£\ 2 ) f{£), where / is the Fourier 
transform of a function / on WL d . It is known that the Levy measure of the process Y has a density 
given by J(x) = j(\x\) where 

/•oo 

j(r) := / (47rt)- d / 2 e- r2/(4t VW*, r > 0. (1.3) 
Jo 

Note that the function r 1— > j(r) is continuous and decreasing on (0, 00). 

We will assume that X is a purely discontinuous rotationally symmetric Levy process with Levy 
exponent v &(^). Because of rotational symmetry, the function ^ depends on |£| only, and by a slight 
abuse of notation we write = ^(|^|). The infinitesimal generator of X is C x := — \P(\/— A). 
We further assume that the Levy measure of X has a density with respect to the Lebesgue measure 
on M. d , which is denoted by Jx{x,y) = Jx{x — y) = jx(\y — %\)- That is, 



with 



e it-(X t -X ) 



e -mm) f or every x eR d &nd££ R d , 



*(l£l)=/ {l-co*(Z-x))Jx(x)dx. (1.4) 
We assume that jx( r ) is continuous on (0,oo) and that there is a constant 7 > 1 such that 

7~ j( r ) — Jx( r ) < 7j( r ) for all r > 0. (1-5) 



This implies that J and Jx are comparable. Clearly (1.5) also implies that 



7-V(l£| 2 )<*(l£l)<70(l£| 2 ) foralUeIR d . (1.6) 

We remark that under the above assumptions, X does not need to be a subordinate Brownian 
motion because jx does not need to be monotone. For example, choose e > such that 2~ 1 j(l) < 
j(l + e) and a continuous function /1 with h(l) = \j{r), h(l + e) = and < h(r) < hj(r) for all 
r > 0. Then jx( r ) '■= j( r ) — h(r) is not monotone and its corresponding Levy process X (through 



Levy exponent (1.4)) is not a subordinate Brownian motion. 
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Under the above setup, X has a continuous transition density p(t, x, y) with respect to the 
Lebesgue measure on M. d (see [9]). Clearly, p(t,x,y) is a function depending only on t and \x — y\, 
and so, by an abuse of notation, we also denote p(t,x,y) by p(t, \x — y\). For every open subset 
D C M. d , we denote by X D the subprocess of X killed upon exiting D. It is known (see )0j ) that X D 
has a transition density PD(t,x,y), with respect to the Lebesgue measure, which is jointly locally 
Holder continuous. Note that pn{t,x,y) is the fundamental solution for C x = — ^(\/— A) in D 
with zero exterior condition and so it can also be called the Dirichlet heat kernel of C x in D. The 
purpose of this paper is to establish sharp two-sided estimates on pjj(t,x,y). The following two 
conditions will be needed for some of the results in this paper when D is unbounded. 

(B) : There exist constants C\ > and C2 E (0, 1] such that 

p(t, u) < Cip(t, C 2 r) for t E (0, 1] and u > r > 0. 

(C) : There exist constants C3 > and C4 E (0, 1] such that 

p(t, r) < C 3 tj(Cir) for t E (0, 1] and r > 0. 



Throughout this paper we will use $ to denote the function 



*(r) = yt~T2\"> r > 0. (1.7) 



ir- 



Note that in particular it follows from (1.2) that 

*( 2r ) = 0(r - 2/4) ^ 0-i( r -2 )/4 =4(& ( r ) for every r>0. (1.8) 

The inverse function of $ will be denoted by the usual notation <3? -1 (r). Here and in the following, 
for a, b £ M, a A b := min{a, 6} and a V 6 := max{a, 6}. 



Remark 1.1 (i) The condition (B) is pretty mild. When X is a rotationally symmetric Levy 
process such that r \-t jx{f) is decreasing, condition (B) holds for all t > (see [251 Propo- 
sition]). In particular it holds for all subordinate Brownian motions with C\ = C% = 1. In 
this special case, we can also see this using the following elementary argument: When X is 
a subordinate Brownian motion, p(t,r) = E[po(<St, r)], where po(t, \x — y\) is the transition 
density of the Brownian motion W . It follows immediately that pit, r) is decreasing in r and 
so (B) holds with C\ = C 2 = 1. 

(ii) Under condition (A), condition (B) is weaker than (C). Under condition (A), we will show 
in this paper that there exists c > 1 such that 

c" 1 ($ _1 (t) d /\tj(r)) <p(t,r) < c(^ 1 (t)) d for t G (0,1] and r > 



(see Proposition |2.2| and Theorem |3 . 7| below) . Thus condition (C) amounts to say that there 
exist constants c > 1 and C4 E (0, 1] such that for i E (0, 1] and r > 0, 

c -1 (* _1 (t) d Atj(r)) <p(t,r) < c($~ 1 (*) d Atj(C 4 r)). (1.9) 



Since j( r ) is a decreasing function in r, (1.9) implies that condition (B) holds with C\ = c 2 
and C2 = C4. 
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(iii) Assume that condition (A) holds. It follows from pU Unj HB] that, for every R > 0, there is a 
constant c = c(T,R,j,<f)) > 1 so that (fijl holds for (t,r) G (0, 1] x (0,i?] with C 4 = 1. (See 



Proposition 2.2 below.) So the assertions in conditions (B) and (C) are always satisfied for 



< r < u < R. 

(iv) By nsansi, 

c- 1 Atj(C 4 -V)) <p(i,r) <c((*- 1 (i)r d Aij(C4r)) for (t,r) G (0, 1] x (0, oo), 

and consequently conditions (B) and (C), hold for a large class of discontinuous processes 
including mixed stable-like processes (with C4 = 1) and relativistic stable-like processes (with 
C 4 = 1, see [HI Theorem 4.1]). □ 



Before stating the main results of this paper, we need first to set up some notations. Let d > 1. 
We denote the Euclidean distance between x and y in U. d by \x — y\ and denote by B(x, r) the open 
ball centered at x G M d with radius r > 0; for any two positive functions / and g, f x g means that 
there is a positive constant c > 1 so that c _1 g < f < eg on their common domain of definition; 
for any open D C M rf and x G D, diam(D) stands for the diameter of D and 5d(x) stands for the 
Euclidean distance between x and D c . 

Definition 1.2 Let < k < 1. VFe say i/iai a open se£ -D is Ac-/a£ «/ i/iere is i?i > suc/i i/iai /or 
all x & D and all r G (0, there is a ball B(A r (x), nr) C D Pi S(x, r). T/ie pair (i?i, k) is called 
the characteristics of the K-fat open set D. 

The following factorization of the Dirichlet heat kernel is the first main result of this paper. 
Recall that C\ and C2 are the constants in condition (B). 

Theorem 1.3 Let X be a purely discontinuous rotationally symmetric Levy process with Levy ex- 



ponent ^ and Levy density Jx satisfying (1.6) and (1.5) respectively, where the complete Bernstein 



function 4> satisfies (A). Suppose that D is a K-fat open set with characteristics (R\,k). 

(i) For every T > 0, there exists c\ = ci(R±, k, 7, T, d,(f>)>0 such that for < t < T, x,y G D, 

PD(t, x, y) > Ci F x (t d > t)F y (T D > t) (^- 1 (t)~ d A U(x, y)) . (1.10) 

(ii) // D is unbounded, we assume in addition that condition (B) holds. For every T > 0, there 

exists C2 = C2(Ci, C2, Ri, k, T, d, 7, <p) > such that for < t <T, x,y G D, 

PD(t, x, y) < c 2 F x (t d > t)F v (r D > t)p(t, C 5 x, C 5 y), (1.11) 

where C 5 = Cf/4. 

(iii) Suppose in addition that D is bounded. Then there exists C3 = c^{diam{D),R\, k,o!, 7,^) > 1 
so that for all (i, x, y) G [3, 00) x D x D, 

% X Wx{t d > l)W y (T D > l)e- tXl < PD (t,x,y) < c 3 F x (t d > 1)P„(td > 1) e~* A \ 

where — Ai < is i/ie largest eigenvalue of the generator of X D . 
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When X is a rotationally symmetric a-stable process in M. d , that is, when = for some 



a S (0, 2), parts (i) and (ii) of Theorem 1.3 are proved in [3]. 

Recall that C3 and C4 are the constants in condition (C). Combining Theorem |1.3[ i)-(ii) and 
Remark |l.l[ ii), we have the following corollary. 

Corollary 1.4 Let X be a purely discontinuous rotationally symmetric Levy process with Levy 



exponent VP and Levy density Jx satisfying (1.6) and (1.5) respectively, where the complete Bern- 



stein function <f> satisfies (A). Suppose that D is a K-fat open set with characteristics (Ri,k). If 
D is unbounded, we assume in addition that condition (C) holds. For every T > 0, there exist 
c\ = c±(Ri, k, T, d, 7, <f>) > and c 2 = 02(63, C4, R\, k,T, d, 7, (j>) > such that for < t < T, 
x,y <E D, 

Ci F x (t d > t)F y (T D > t) ($>-\t)- d A tj(\x - 
<p D (t,x,y) < c 2 F x {t d > t)F y (T D > t) ($- 1 (ty d Atj(C e \x - y|)) , 

where C 6 = C|/4. 

The second main result of this paper is on explicit sharp Dirichlet heat kernel estimates for 
subordinate Brownian motions in C ' open sets. So in the remainder of this section, we assume 
that X = Y, a subordinate Brownian motion with Levy exponent Vl/(£) = </>(|£| 2 ). 

Recall that an open set D in M. d (when d > 2) is said to be a (uniform) C 1,1 open set if 
there exist a localization radius R2 > and a constant A > such that for every z £ dD, there 
exist a C^-function tp = ip z : E d_1 -> R satisfying V(0) = 0, W(0) = (0,...,0), ||V?/>||oo < A, 
^^(x) — VV'(^)! < A|x — z|, and an orthonormal coordinate system CS Z with its origin at z such 
that 

B(z, R 2 )C\D = {y= (y, y d ) in CS Z : \y\ < R 2 ,y d > rp(y)}. 

The pair (R2, A) is called the characteristics of the C > open set D. Note that a C 1 ' 1 open set D 
with characteristics (i?2j A) can be unbounded and disconnected; the distance between two distinct 
components of D is at least R 2 . By a C 1 ' 1 open set in M we mean an open set which can be written 
as the union of disjoint intervals so that the minimum of the lengths of all these intervals is positive 
and the minimum of the distances between these intervals is positive. 

Here is the second main result of this paper, which gives an affirmative answer to the Conjecture 
posed in |16j . In view of Remark |1.1[ it extends the main results of jl 1|, I13 |. [T4^| . Recall that condition 
(B) holds with C\ = C 2 = 1 for any subordinate Brownian motion. 

Theorem 1.5 Suppose that X is a subordinate Brownian motion with Levy exponent = 
<Ml£| 2 ) with (j) being a complete Bernstein function satisfying condition (A). Let D be a C 1,1 open 
subset of~R d with characteristics (R 2 ,A). 

(i) For every T > 0, there exists c\ = c\{R 2 , A, T, d,4>) > such that for all (t, x, y) S (0, T]xDxD, 
m[t , x , y) > ci (l A «) V2 (l A i^Mf (.-.<*)- A «<,,,)) . 
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(ii) For every T > 0, there exists C2 = C2(i?2, A, T, d, 4>) > such that for all (t,x,y) G (0,T] x 
DxD, 

PD(t , x, „ < C2 fl A 5%M1) "Vl A ™) 1/2 P (, |, - ,|/4). 



(iii) Suppose in addition that D is bounded. For every T > 0, there exists C3 > 1 depending only 
on diam(D), A, R2, A, d, and T so that for all (t, x, y) G [T, 00) x D x D, 



cj 1 e- Xlt y/*(8D(x))y/$(5 D (y)) < p D (t,x,y) < c 3 e~ Al V$(M*))\/$(Mi/))> 
where —Xi<0is the largest eigenvalue of the generator of X D . 



Recall that C3 and C4 are the constants in condition (C). Combining Theorem |l,5[ i)-(ii) and 
Remark |l.l[ ii), we have the following corollary. 

Corollary 1.6 Suppose that X is a subordinate Brownian motion with Levy exponent = 
0(|^| 2 ) with (f> being a complete Bernstein function satisfying condition (A). Let D be a C open 
subset ofW 1 with characteristics (i?2,A). If D is unbounded, we assume in addition that condition 
(C) holds. For every T > 0, there exist c\ = ci(i?2> A, T, d, 0) > and C2 = 02(6*3, C4, -R2, A, T, d, 0) > 
smc/i that for < i < T, x, y G D, 

ci ^ A " ^ A " 2 ( ,_ I(tr , A t _ s|)) 

< M*. ») < <* (l A (l A V2 (.-'(*)-" A «C 4 |* - 9 |/4)) . 

When X is a rotationally symmetric a-stable process in M. d , Theorem 1.5 is first established in 
[TTj . Sharp two-sided Dirichlet heat kernel estimates in C 1 ' open sets are subsequently established 
in [12l HH [131 E] f° r censored stable processes, relativistic stable processes, mixed stable processes, 
and mixed Brownian motion and stable processes, respectively. By integrating the two-sided heat 
kernel estimates in Theorem 1 1 . 5| with respect to t, we obtain the two-sided estimates on the Green 



function Gn(x,y) := J °° p£>(t,x,y)dt (see Theorem 7.3 below), which extend [201 Theorem 1.1]. 



Condition (A) is a very weak condition on the behavior of (f> near infinity. Using the tables 
at the end of [23], one can come up plenty of explicit examples of complete Bernstein functions 
satisfying condition (A). Here are a few of them. 

(1) (f>(\) = A a / 2 , a G (0,2] (symmetric a-stable process); 

(2) 0(A) = (A + m 2 / a ) a / 2 — m, a G (0, 2) and m > (relativistic a-stable process); 

(3) 4>(X) = \ a l 2 + A' 3 / 2 , < P < a < 2 (mixed symmetric a- and /3-stable processes); 

(4) <j>(\) = A a / 2 (log(l + \))P, a G (0, 2), p G [-a/2, (2 - a)/2]. 

Now we give a way of constructing less explicit complete Bernstein functions that have very 
general asymptotic behavior at infinity. Suppose that a G (0, 2) and £ is a positive function on 
(0, 00) which is slowly varying at infinity. We further assume that t -»■ t a l 2 l{t) is a right continuous 
increasing function with lim^o t a / 2 £(t) = (so Jq°°(1 + t)~ 2 t a / 2 £(t)dt < 00). Then the function 



/•oo 

/(A) := / (A + t)- 2 t a / 2 £{t)dt 
Jo 
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is a Stieltjes function, and so the function 



0(A) := J- = n°°(X + trH^£(t)dt\ 



is a complete Bernstein function (see |23l Theorem 7.3]). It follows from |26|, Lemma 6.2] that 
0(A) x A Q / 2 ^(A) when A > 2. 

The rest of the paper is organized as follows. Section [2] recalls and collects some preliminary 
results that will be used in the sequel, including on-diagonal heat kernel estimates and the boundary 
Harnack principle. Section [3] presents the interior lower bound heat kernel estimates, including an 
explicit lower bound estimate for symmetric Levy processes on M. . The proof of the short time 
factorization result for prj(t,x,y) (that is, Theorem |1.3[ i) and (ii)) is given in Section|4j while the 
proof of Theorem |1.5| (i) and (ii) is given in Section [5| The large time heat kernel estimates are 
proved in Section [6j The Green function estimates for subordinate Brownian motions in bounded 
C 1 ' 1 open sets are derived in Section[7]from the two-sided Dirichlet heat kernel estimates in Theorem 



1.5. The derivation, however, requires quite some effort. 

Throughout this paper, d > 1 and the constants Ro, Ri, R2, A, k, Si, 62, C\, C2, C3, C4, C5 
and Cq will be fixed. We use ci,C2,--- to denote generic constants, whose exact values are not 
important and can change from one appearance to another. The labeling of the constants ci, C2, • • • 
starts anew in the statement of each result. The dependence of the constant c on the dimension 
d will not be mentioned explicitly. We will use ":=" to denote a definition, which is read as "is 
defined to be". We will use d to denote a cemetery point and for every function /, we extend its 
definition to d by setting f{d) = 0. We will use dx to denote the Lebesgue measure in M d . For a 
Borel set A C M. d , we also use \A\ to denote its Lebesgue measure. 

2 Preliminary 

In the first part of this section, we assume that X is a purely discontinuous rotationally symmetric 



Levy process with Levy exponent ^ and Levy density Jx satisfying (1.6) and (1.5) respectively, 
where the complete Bernstein function satisfies (A). 

A function u : M rf 1— > [0, 00) is said to be harmonic in an open set D C M rf with respect to X if 
for every open set B whose closure is a compact subset of D, 

u(x) = E x [u(X TB )\ for every x G B. (2.1) 

A function u : K 1— > [0, 00) is said to be regular harmonic in an open set D C M. d with respect to 
X if 

u(x) = K x [u(X TD )] for every x G D. 

Clearly, a regular harmonic function in D is harmonic in D. 

Very recently the following form of the boundary Harnack principle is established in |21j . 



Theorem 2.1 ([21, Theorem l.l(i)]) There exists c = c((p,j) > such that for any z$ G 
any open set D C M d , any r G (0, 1) and any nonnegative functions u, v in ~R d which are regular 
harmonic in D D B(zq, r) with respect to X and vanish in D c n B(zQ,r), we have 

^rl <c^H for allx,y G Dn%r/2). 
v(x) v{y) 
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Recall also from |21j that j enjoys the following properties: for every F > 0, 



j(r) 



fd <])(j~) 



for r G (0,F], 



(2.2) 



and there exists c > such that 



j(r) < cj(r + 1) 



for r > 1. 



(2.3) 

It follows from (1.8) and (2.2) that both the function <3? defined by (1.7) and the function j satisfy 



a doubling property; that is, for every constant R > 0, there is a constant c > 1 so that 

$(2r) < 4 3>(r) and j(r) < cj(2r) for every r G (0,F]. 

Moreover, under condition (A), for any given R > 0, j(r) satisfies all the conditions in 
r G (0, R). Thus we have the following two-sided estimates for p(t, x, y) from [91 HOj : 

Proposition 2.2 For any T > 0, there exists c\ = ci(T,R,j,(f>) > such that 

p(t, x, y) < ci ($-\t))- d for (t, x, y) G (0, T] x R d x M d . 

For any T,R> 0, there exists C2 = C2(T, F, 7, 0) > 1 snc/i that for all (t,x,y) G [0, T] x 1 
wrat/i |x — y| < F, 

C 2 1 ((d)- 1 (t))- d A t J(x, y)) < p(t, x, y) < c 2 (($- 1 (i)) _d A i J(x, y)) . 



(2.4) 
181 for 



(2.5) 



i d x R d 



(2.6) 



Proof. (2.5) is given in the first display on page 1073 of [Sj- By [9j Theorem 2.4] and (1.5), there 
exist F* > and F* such that for all (t, x, y) G [0, T«] xR d xR d with \x-y\ < F*, ((2~6) holds. Now 



we assume F > F*. We construct Z from X by removing jumps of size larger than F via Meyer's 
construction (see [22]). Let pz(t,x,y) be the transition density of Z. By [U Lemma 3.6] and [2J 
Lemma 3.1(c)] we have for every t > and x, y G 



3 _ *|| JiilloO 



where 



p z (t,x,y) < p(t,x,y) < p z (t,x,y) + t\\J R \\ c 



Jn(x,y) := Jx(x,y)l{\ x - y \ >R } and Jr(x) := J R (x,y)dy. 



(2.7) 
(2.8) 



Applying |10l Theorem 1.4] and its proof to (2.7), we get 

c\e~ 



^-\t))- d MJ x (x,y)) <p(t,x,y) < c 2 ((^ 1 (t))- d A t J x (x, y)) + t\\ J^. (2.9) 



(2.6) now follows from 



□ 



The function Jx(x,y) gives rise to a Levy system for A, which describes the jumps of the 



process A: for any non- negative measurable function / on 



x 

y G M. d and stopping time T (with respect to the filtration of A) 

f(8,Xs-,X a ) =E X f 

s<T 



x R d with f(s, y, y) = for all 



Ex 







f(s,X s ,y)J x (X s ,y)dy ds 



(2.10) 
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(See, for example, [HI Proof of Lemma 4.7] and [THl Appendix A].) 

When = <?K|£| 2 ) an d <f> is a complete Bernstein function satisfying condition (A), the 

following boundary Harnack principle on a C ' open subset with explicit decay rate is established 
in HI] (see also [20]). 

Theorem 2.3 ([19, Theorem 1.5]) Suppose that X is a rotationally symmetric Levy process 
with Levy exponent ^{Cj = <?K|£| 2 ) with 4> being a complete Bernstein function satisfying condition 
(A). Assume that D is a (possibly unbounded) C 1 ' 1 open set in R d with characteristics (1?2,A). 
Then there exists c = c(l?2, A, <p, d) > such that for r G (0,(122 A l)/4], Q G dD and any 
nonnegative function u in M. d that is harmonic in D H B(Q,r) with respect to X and vanishes 
continuously on D c n B(Q,r), we have 



u(x) 



< c 



u(y) 



for every x,y G D n B(Q, r/2). 



(2.11) 



*(fe(x))V2 " $(5 D (y))V2 

3 Interior lower bound estimate 

In this section, we assume that A is a purely discontinuous rotationally symmetric Levy process with 
Levy exponent ^ and Levy density Jx satisfying (1.6) and (1.5) respectively, where the complete 
Bernstein function <p satisfies (A). We will give some preliminary lower bounds on p B (t,x,y) and 
p(t,x,y). We first recall the following from (9j[18]. Recall that J B was defined in (2.8). 

Lemma 3.1 For any positive constant M, there exist e = e(M, 7,0) > and c = c(M,7, (f>) > 1 
such that for all r G (0, M] and z G M d , 



i P -\\Ji\l 



and 



{tB{z,t) > £ $( r )) > 2 e 



'^(r) < E z [t bM ] < c*(r). 



(3.1) 
(3.2) 



Proof. (3.1) follows directly from |9, Lemma 2.5]. We then have 

IE, [r B (z,r)] > e<S>(r)F z (T BM > e*(r)) > 2- 1 e~^°°e& 



r). 



On the other hand, by the Levy system for X in (|2.10|), (|1.5|), (|2.2|) and (|2.4|), we have 
1 > ¥ z (X TbM G 5(z,2r) c ) = E 2 



T B(z,r) 



B(z,2r) c 



J x (X s ,y)dyds 



> ciE 2 



T B(z,r) 



'B(je,3r)\B(j;,2r) 
which yields E 2 [t B ( 2jJ ,)] < c<&(r). 



J(X s ,y)dyds 



" C2 E * L r S(*.r)J 



□ 



Lemma 3.2 For any positive constants a and R, there exists c = c(a, 12, 7, </>) > suc/i i/jai /or aW 
2 G M d and r G (0,12], 

/9 ,M r B(*,r) > a ®(r)) > c 
y£B{z,r/2) 
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Proof. By (3.1) and (2.4), there exists e± = £\{R, 7,0) > such that for all r G (0, R], 

inf F z (T B(z , r/2) > e x *(r)) > 2-V" ^H°° . 



Thus it suffices to prove the lemma for a > e\. Applying the parabolic Harnack inequality 
Theorem 5.2] at most 2 + [a/e{\ times, we conclude that there exists c\ = ci(a, R, 7, (ft) > such 
that for every w,y G B(z,r/2), 



Thus 



(TBfer) > a$(r)) = / p B{z . r) (a<S>(r),y,w)dw 

JB(z,r) 



l(z,r) 

> 



PB(z,r)(a®(r),y,w)dw 

B(z,r/2) 



> ci p B , z>r / 2 )(£x<$>(r),z,w)dw 

JB(z,r/2) 

= ciP z (r B( , ir/2) > £l $(r)) > 2- 1 e-ll^ll°°. 
This proves the lemma. □ 

For the next four results, D is an arbitrary nonempty open set and we use the convention that 
<$£>(■) = 00 when D = R d . 

Proposition 3.3 Let T > and a > be constants. There exists c = c(T, a, 7, (/>) > suc/i i/iai 

PD(t,x,y) > c($- 1 (t))- d (3.3) 
/or every (t,x, y) € (0, T] x D x D wt/i 5d(x) A <5o(y) > a<I> _1 (t) > 4\x - y\. 

Proof. We fix (t,x,y) £ (0,T\ x D x D satisfying S D (x) A S D (y) > a®' 1 ^) > i\x - y\. Note that 
\x-y\< a<5>- l {t)/A < a<5>- l (T)/A and that 

B{x,a$- l {t)/4) C B(y,a$ _1 (t)/2) C B(y, 2a^ l {t)/3) C L>. 

So by the parabolic Harnack inequality (9j Theorem 5.2], there exists c\ = ci(T, 7, cp) > such that 

ciPD(t/2,x,w) < pn-(t,x,y) for every w G B(x, a<J> -1 (r)/4). 



This together with Lemma 3.2 and (|2.4[) yields that 



p D (t,x,y) > — ^fTTs / -j \ I / p D (t/2,x,w)dw 

\B(x,a<£> l (t)/A)\ J B (x,a^-Ht)/4) 

> c 2 (<S>-Ht))~ d / p B(x ^-i {t)/4) (t/2,x,w)dw 

JS(x,a'I>- 1 (t)/4) 

= C 2 (*-Ht))- d Vx(T B ( x ,a*-Ht)/4) >t/2) > C 3 (*-\t))- d , 

where Cj = Cj(T, a, 7, 0) > for i = 2,3. □ 
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Lemma 3.4 Let T > and a > be constants. There exists c = c(a, T, 7, <p) > so that 

F x (X? G o*" 1 ^)^)) > c($- 1 (t)) d tJ(x,y) 
/or every y) G (0, T] x D x D mt/i <5d(x) A 5d{u) > a<5 _1 (i) and a<I> _1 (t) < 4|x - y|. 



Proof. It follows from Lemma 3.2 that, starting at z G -B(y, a<I> _1 (i)/4), with probability at least 
ci = ci(a, T, 7,0) > the process X does not move more than a<I> _1 (i)/6 by time t. Thus, it 
suffices to show that there exists C2 = 02(0, T, 7, 4>) > such that for (t, re, !/)e(0,T]xflxD with 
A 5 D (y) > a$- 1 (i) and a* -1 ^) < 4\x - y|, 



? x (X D hits the ball £(y, a$ -1 (*)/ 4 ) by time t) > c 2 ($~ 1 (t)) d tJ(2;,y). 



(3.4) 



Let := B(x, a<I> 1 (t)/9), £?* := B(y, a® 1 (t)/9) and r* := r B t. It follows from Lemma 



and (2.4) that there exists C3 = 03(0, T, 7, <fi) > such that 

Ex [t A t*] > t Px (r* > i) > c 3 t for all t G (0, T]. 



3.2 



(3.5) 



Since B x n -B* = 0, by the Levy system of X and (1.5), 

Px (X D hits the ball B(y, a<&~ l (t)/4) by time t) 



> Px (X tAr t G B(y, o*- 1 ^)/^) > E 2 



tArf 



1/ 



Jx(X s ,u)duds 



tArl 



y 



J(X S , u)duds 



(3.6) 



We consider two cases separately. 

C aS" 1 ^). Si 

|X S - u| < I X s - rc| + I re - y\ + |y - u\ < 2\x - y\. 



(i) Suppose |rc — y| < a*!* 1 (T). Since \x — y\ > a& (t)/4, we have for s < r* and it G By, 



Thus by (3.5) and (3.6), 



F x (X D hits the ball B(y, a$ _1 (i))/4 by time t) 
> C4Mx[tAri] \B y \j(2\x-y\) > c 5 (t)) d tj (2\x - y|) 



for some positive constant C5 = c^(a, T, 7, 0) > 0. Therefore, in view of (2.4), the assertion of the 
lemma holds when |rr — y| < a3> _1 (T). 

(ii) Suppose \x — y| > a$ _1 (T). In this case, for s < r * and u G B y , 

\X S - u\ < \X S -x\ + \x-y\ + \y-u\<\x-y\+ a$ _1 (*)/ 4 < \ x ~ v\ + a$ _1 (T)/4. 



Thus from (3.6) and then (3.5), 



F x (X D hits the ball B(y, a^ l {t)/A) by time t) 
> cqE x [t A r*] / j (|x-y| +a$- 1 (T)/4) 



f/ 



> 07*15*1^(1^-2/1 + a*- 1 (r)/4) 

> CBt^-Ht^dx-yl+a*- 1 ^) 
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for some constants q = Cj(a, T, 7, (f) > 0, i = 6,7,8. Thus we conclude from (2.3) and (2.4) that 
the assertion of the lemma holds for \x — y\ > a<I>~ 1 (T) as well. □ 



Proposition 3.5 Let T and a be positive constants. There exists c = c(T, a, 7,0) > such that 

p D (t,x,y) > ctJ(x,y) 
for every (t,x, y) G (0, T] x D x D with 5d{x) A 5d{u) > a$ _1 (i) and a^ 1 ^) < A\x - y\. 



Proof. By the semigroup property, Proposition 3.3 Lemma 3.4 and (2.4), there exist positive 



constants c\ = ci(T, a, 7, 4>) and C2 = C2(T, a, 7, 0) such that 



PD(t,x,y) 



PD(t/2, x, z)p D (t/2, z, y)dz 



D 

> / PD(t/2,x,z)p D (t/2,z,y)dz 

JB(y,a^{t/2)/2) 



> Cl {^-\t/2))- d ¥ x \Xf /2 G B(y, a^-\t/2)/2) 

> c 2 tJ(x,y). 



□ 



Combining Propositions 3.3 and 3.5 , we have the following preliminary lower bound for pu{t, x,y). 



Proposition 3.6 Let T and a be positive constants. There exists c = c(T, a, 7,0) > such that 

PD(t,x,y) > c((^-\t))- d AtJ(x,y)) 
for every (t,x,y) G (0,T] x D x D with 5 D (x) A5 D (y) > a$ _1 (t). 

id 



In particular, Proposition 3.6 



with D 



gives 



Theorem 3.7 For any constant T > 0, there exists c = c(T, 7, (f>) > such that for all (t,x,y) G 
[0, T] x R d x M d , 

p(t,x,y) >c ((<5>- 1 {t))- d AtJ(x,yf 



4 Factorization of Dirichlet heat kernel 

In this section, we continue assuming that X is a purely discontinuous rotationally symmetric Levy 



process with Levy exponent ^ and Levy density Jx satisfying (1.6) and (1.5) respectively, where 



the complete Bernstein function (ft satisfies the growth condition (A). 

Throughout this section, T > is a fixed constant and D is a fixed K-fat open set with char- 



acteristics (Rx,k). Recall that A r {x) G D is defined in Definition 1.2 For (t, x) G (0,T] x D, set 
r = r(t) = S-^fli/^-^T) < Ri and define 



U(x,t) := DH B(x, \x - A r (x) \ + Kr/3), V(x,t) := DnB(x, \x - A r (x)\ + kx). (4.1) 
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Figure 1: U(x,t) and V(x,t) 



Let 

A' r (x) £ D be a point such that B(A' r (x), nr/S) C B(A r (x),nr) \ U(x,t). (4.2) 

Note that B(A r (x), nr/3) cU(x,t) and B(A' r (x), rer/6) CB(^.(4/tr/3) C t) \ t). See 
Figure [T} 

Lemma 4.1 For ewery T > and M > 1, we /iaue t/iat, for (t,x) 6 (0, T] x D, 

P x (r D > t/Af) x P^Ty^ > Aft) x F x (T V{x<t) > t/M) x F x {t d > Aft) 

x F x (X Tu(x t) eD)^ rXlr^tj], (4.3) 



where U(x,t) and V(x,t) are the sets defined in (4.1) and t/ie (implicit) comparison constants in 



(4.3) depend only on d, M,T, R\, k, 7 and 



Proof. Without loss of generality we assume f?i < 1. Fix (t,x) £ (0, T] x D, and set r = r(t) = 
$ -1 (t).Ri/$ -1 (T) < i?i < 1. Recall that U(x,t) and F(x,i) are the sets defined in fl3~I] ). Observe 
that 

Pa.CTVC.t) > Aft) < P^Ty^) > t/M) A P x (t d > Aft) < P a .(r D > t/Af). (4.4) 
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Note that by (O), (272b and (12.101), we have 



T U(x,t) 



v (X(r u{X)t) )eB(A' r (x),Kr/6)) =E X I I J x (X t ,y)dtdy 

JO JB(A' r (x),Kr/6) 

r d j(r)E x [T U{x , t) ] x ^r^M^ru^] x t" 1 E x [t u{x ^}. 



(4.5) 



If | a; — A r (x)| < Kr/2, then -B(x, nr/3) C {/ (x, t) C V(x, t) and so 

1 > W x (td > t/M) > F x (t d > Mt) > F x (T V(Xtt) > Mt) > F x (r B{XjKr/3) > Mt), 
which is greater than or equal to a positive constant depending only on 0, T, R\, 7, d, k, M by 



Lemma 3.2 Thus we have, in view of (4.4) 



\{t d > t/M) x F x (T V(Xtt) > Mt) x F x {t v(xA > t/M) x F x (t d > Mt). 



Moreover, when \x — A r (x)\ < Kr/2, we have B(x, nr/3) C U(x, t) C B(x, r) and so by Lemma 3.1 



and (2.4), 



cit < c 2 $(Kr/3) < E x [T B{xKr/3) ] < E x [t u{x ^} < E x [T B{Xtr) ] < c 3 $(r) < c A t. 



Combining the last two displays with (4.5) and the fact that 



1 > MXr v{x , t) eD)> F x {X Tv{xt) G B(A' r (x),Kr/6)), 

we arrive at the assertion of the lemma when \x — A r (x)\ < Kr/2. 
Now we assume that \x — A r (x)\ > Kr/2. We note that 

F x {t d > t/M) < F x (T U{x>t) > t/M)+F x {X Tu{xt) G D) . 

For r G (0, Ri], by Theorem |2.1| we have 

\{X Tu(xt) eB(A' r (x),Kr/6)) 



(4.6) 



V(x)K (l , () eB(A' r (x), K r/6)) 



< c 5 - 



G 5(^;(x),Kr/6)) 



J 4 r (x) (-^rB(A r (a:),K)-/3) 



G £(A;(x),Kr/6)) 



Note that when (w, y) G B(A r (x), nr/3) x 5(^4' r (x), «x/6), — j/| < 2w and so — y|) > j(2nr) 



Thus, by ([1_5|), ( |2.10p and Lemma |3.1[ 

€B(A' r (x),Kr/6)) =E Ar{x) [ 

JBIAL 



A r (x) ip^T B ( Ar ^ x) Kr /^ 



Jx(X t , z)dtdz 



B(A' r (x),Kr/6) JO 

> c Q r d j(r)E Mx) [T B ( Mx)jKr/3) ] > c 7 $(r)~ 1 Eo[r B(oW3) ] > cg^r)" 1 $(/cr/3) > eg, 
where the last inequality is due to ( |2.4[ ). Therefore, 

M X tu^£D) < c 10 F x {X Tu{xt) €B(A' r (x),Kr/Q)). 



Since F x {r u{Xt t) > t/M) < Mt~ 1 E a! [ru(x,t)] , w e have by (4.6), (|47j) and j45] ), 

> t/M) < cut-X^^t)]. 



(4.7) 



(4.8) 
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On the other hand, by the strong Markov property of X, 

M X T U(x<t) £B(A' r (x),Kr/6)) 
= E, [Px^ (rB(X Tu(xtyK r/e) > Mt) ;X Tu ^ t) G B(A' r (x), K r/6) 
+ E - [ P ^ ( .,i) ( r B(^ (ljt)lW /6) < Mt) : X Tt7(M) G B(^(i),«r/6) 

+P (t b((W 6) < Mt) P x (x T[/(:M) G B(4(x), kt/6)) 
= P, (ry^t) > Mt, X T[/(M) G B(A' r (x), «r/6)) 

+ (1 - P (r B(0 ,« r/6 ) > Mt)) P x (-X^, G S(4( S ),Kr/6)) 
< P, (r^ lt) > Mt) + (1 - P (r B((W6) > Mt)) P x G w/6) 

Thus 

Po(r B((W 6) > Mt)P x (X^ ( , it) G B(A' r (x),Kr/6)) < F x (r v{x , t) > Mt) . 
Since Pq (tb(p K r/6) > Mt) > C12 by (2.4) and Lemma 



(4.9) 



3.2 



combining (|4.5[) with (4.9), we get 



r^tTfrOM)] < ci 3 P,(X Tc/(;Cit) G B(A' r (x),Kr/6)) < c 1A P x {T V{x ,t) > Mt). 



This together with (4.8) and (4.4) completes the proof of this lemma. 



□ 



Lemma 4.2 Suppose that U\, U3, E are open subsets ofM. d with U\, Us C E and dist(C/i, C/3) > 0. 
Let U 2 := E\ (U\ UU 3 ). If x G U x and y G U 3 , then for all t > 0, 

PE(t,x,y) <P X (X TU1 G U 2 ) f sup p E (s, z,y)\ + (t A E x [rc/J) ( sup Jx(u,*)) (4.10) 



and 



PE(t,x,y) > tP x (Tu 1 >t)Py(ru a >t) inf J* («,*). 

U&J1.Z&J3 



(4.11) 



Proof. For (4.10), see the proof of |15| Lemma 3.4]. For (4.11), see the proof of |16} Lemma 3.3]. 
□ 



Proof of Theorem 1.3| (ii). It follows from Lemma 4.1 and the semigroup property that it suffices 
to establish the assertion for T < 1. So in the remainder of this proof, we assume that T < 1. Fix 



t G (0,2*1 and set r = $- 1 (t)iii/$- 1 (T). By dO), fOl), D) and Theorem 



3.7 



p(t/2, x, y) x (^(t))-* if \x-y\< Sr. 



By the semigroup property, (2.4), (2.5) and Lemma |4.l[ when \x — y\ < 8r, 

p D (t/2,x,y) = / p D (t/4:,x,z)p D (t/4,z,y)dz 



D 



< swpp(t/4,z,y)P x (T D >t/4) 



z& 



< ciC*" 1 ^))-^^ > t). 



(4.12) 



(4.13) 
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Thus by (14121) and CTl, 



p D (t/2, x, y) < c 2 P x (t d > t)p(t/2, x, y) if \x - y\ < 8r . 



(4.14) 



Now we assume \x — y\ > 8r. Let D\ := U(x,t) be the set defined in (4.1), D 3 := {z G D : 
1 2 — x| > |x — y\/2} and 

D 2 :=D\ [Px UD 3 ) = {z£D\ U(x, t) : \z - x\ < \x - y\/2}. 

Then by condition (B) we have 

sup p(s,z,y)< sup p(s,\z- y\) < C\ sup p(s,C 2 \x- y\/2). 

s<t/2,z&D 2 s<t/2, \z-y\>\x-y\/2 s<t/2 

Extend the definition of p(t, r) by setting p(t, r) = for t < and r > 0. For each fixed x, y G M rf 
and t > with |x — y\ > 8r, one can easily check, in view of (1.9), that (s, w) \-t p(s, C 2 \w — y\/2) 
is a parabolic function in (— 00, T] x B(x,2r). So by the parabolic Harnack inequality from [9j 
Theorem 5.2], there is a constant c 3 = C3(d, 7, 0) > 1 so that for every t G (0,T], 



sup p(s, C 2 \x - 2/|/2) < c 3 p(t/2, C 2 |x - y|/2). 

s<t/2 



Hence we have 



sup p(s,z,y) < c 3 C lP (t/2,C 2 \x-y\/2) < c 3 Cfp(t/2, C%\x-y\/4), 

s<t/2, z&D 2 



(4.15) 



where in the last inequality, we applied condition (B) again. If u E D±, then |x — u\ < \x — A r (x)\ + 
nr/3. Thus if u G D%, z G D 3 , then 

|u — z\ > |z — x| — |x — n| > |z — x| — |x — A r (x)| — nr/3 > \z — x\ — r > -\z — x\ > — \x — y\. 



Thus by (2.2), (1.5), (2.4) and Theorem 3.7 



t sup Jx(u,z) < sup j(\u - z\) < jtj(\x - 2/1/4) < jtj(C 2 \x - y\/4) 

u&D^^&Di \u-z\>\x-y\/i 

< c iP (t/2,Cl\x-y\/A). 



Let C5 := C2/4. It follows then from (4.15), Lemmas 4.1 and 4.2 that for |x — y\ > 8r, 



PD {t/2,x,y) < c 5 p(t/2,C 5 x,C 5 y) (F x (X Tu(x t) G D) + i^E* [7^ 
< c 6 Px(TD>t)p(t/2,C , 5 s J C 5 y). 



(4.16) 



Hence by condition (B), ( |4.14[ ), (4.16), symmetry, the semigroup property and Lemma |4.1| we 
conclude that for all (i, x, y) G (0, T] x D x L>, 



p D (t,x,y) 



D 



PDit/2, x, z)p D (t/2, z, y)dz 



< 4P x {t d > t)F y (r D > t) I p(t/2, C 5 x, C 5 z)p(t/2, C 5 z, C 5 y)dz 



c 7 F x {t d > t)¥ y (T D > t) / p(t/2, C 5 x, z)p(t/2, z, C 5 y)dz 
c 7 F x (r D > t)¥ y (r D > t)p(t, C 5 x, C 5 y) . 
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□ 



Proof of Theorem [Dsfr). Fix (t,x,y) G (0,T] x D x D and set r = $ _:L (t)i2i/$- 1 (T). Let 
U(x,t) be the set defined in (4.1), A r (x) G D and >Jy(x) 6 D be the points defined in Definition 



1.2 and \A.2\ respectively. Note that 
PD(t,x,y) = 



PD(t/3, x, u)pu(t/3, u, v)pi)(t/3, u, y)dudv 



DxD 

> inf p£f(t/3,u,v) 

(u,v)£B(A' r (x),Kr/6)xB(A' r (y),Kr/6) 



B(A' r (x),Kr/6) 



poit/3, x, u)du 



Jb< 



p D (t/3,v,y)dv. 

lB(A' r {y), K r/6) 
For (u,v) G B(A' r (x),Kr/6) x B(A' r (y), «r/6), 

inf p D (t/3,u,v) 

(u,v)eB(A' r (x),Kr/6)xB(A' r (y),Kr/6) 



(4.17) 



> ci inf 

(u,v)£B(A' r (x),Kr/6)xB(A' r (y),Kr/6) 

> c 2 (tJ(x,y)A(^- 1 (t))- d ), 



^-\t/3))- d A(t/3)J(u,v] 



(4.18) 



where in the first inequality we used Proposition 3.6, in the second inequality we used (2.2) and 
(2.3) respectively in the cases \x — y\ > nr and \x — y\ < nr. 



On the other hand, for u G B(A' r (x), Kr/6), by Lemma 4.2 with U\ = U(x,t) and U 3 = 
B(A' r (x), nr/4), we have 

PD (t/3,x,u) > ff,(T %() > t/3)F u (r Uz > t/3) inf J x {w,z) 

weu(x,t), z&u 3 

> t/3) inf j(\w — z\). 
weu(x,t),zeu 3 

Since 

j(\w-z\) >j(r) =j($- 1 (t)R 1 /<S>-\T)) for (w,z) G U(x,t) x U 3 , 
we have by Lemma 3.2, (2.2) and ( |2.4[ ), 

p D {t/3,x,u) > c 3 t¥ x (r u(X:t) > t/3) $-1 * = ^-\t)- d F x {r u{X;t) > t/3). 

It follows from condition (A) that there exists C4 > 1 such that $ _1 (6t) < ab 1 ^ 2 ^^ 1 (t) for every 
t < T and b G (0,1]. Thus with a := (k/(3c 4 )) 2<52 we have that 3~ 1 K$~ 1 (t) > $ _1 (ai) for every 
t <T, and hence V(x,at) C U(x,t). Thus Lemma 4.1 implies that 

Vx(Tu(x,t) > t/3) > Wx(T V (x,at) > t/ 3 ) ^ c 5 F^.(-rz? > t). 



Consequently, by (2.4), 



C6 



p D (t/3,x,u)du > , lu . d 

B{A' r (x), K r/6) * (t) a 



IPx(td > t)|5(4(x),Kr/6)| > ctP^td > t). 



Similarly, using the symmetry we also have f B , A , i y \ Kr / 6 \ po(t/3, v, y)dv > csF v (td > t). Hence we 
have by §4~T7\\ and (|IT8|), 



PD {t,x,y) > c 9 F x (t d > t)F y {r D > t) ((*- 1 (t))- rf AtJ(x,yj) . 



□ 
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5 Small time Dirichlet heat kernel estimates for subordinate Brow- 
nian motions in C 1 1 open set 

In this section we assume that D is a C 1 ' 1 open set in M. d with characteristics (R 2 ,X) and that 
X is a subordinate Brownian motion with Levy exponent = 0(|^| 2 ), where <p is a complete 

Bernstein function satisfying condition (A). 



Proof of Theorem 1.5| (ii). Without loss of generality we assume R 2 < 1. In view of Theorem 
1.3[ ii), it suffices to show that 

1/2 

P* {td >t)<c(lA _^MA for (t } x ) e (0, T] x D. (5.1) 



Fix (t,x) G (0,T] x D and set r = r(t) = $- 1 (t)R 2 /$- 1 (T) < R 2 < 1. By Theorem [3^1 we only 
need to show the theorem for 5d(x) < r/16. Take xq G 9-D such that <Sd(x) = |x — xo|. 

Let U\ := B(xo,r/8) fl .D and n(xo) be the unit inward normal of dD at the point xq. Put 

r 

lb 



Note that 8d(xi) = r/16. Applying the boundary Harnack principle (Theorem 2.3) and (2.4) we 
get 



» x (X TUi eD\Ui)< c^ Xl (X TUi eD\Ui)< 



1 H5 D (x)) 

Hs D (xi)) 



< C ^(X mi eB\uJ^< c J*^ 



(5.2) 



Take x 2 G R d so that B(x 2 ,r) C B(x ,4r)\ B(x ,r). Then, by (2.10), (2.2) and (2.4), we have 



F x (X TUi G B(x 2 ,r)) = E x 



J(\X s -y\)dyds 

JB( X2 ,r) 

> c 3 \B(0, $- 1 (t))|j(5r)E x [r c/l ] 

> c 4 $($- 1 (t))- 1 E :c [r c/l ] = c^-XK]. 



Now by the same argument as that in (5.2), we get, 



E*fa/i] < cJ 1 tP a; (X n;i G B(ar 2 ,r)) < cgtP^ (JC^ G £?(x 2 , r)) J < c 5 $(<5 D (x)). (5.3) 



Thus, by (5.2) and (5.3), we have 



(r D > t) < F x (t Vi >t) + P, (X TUi eD\U 1 



< -E x [rc/J + P x I X Tyi G D \ C/i j < c 6 



$(<5d(x)) 



Since P^ (td > £) < 1, (5.1) follows immediately. 



(5.4) 

□ 
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Let 5qd(x) be the Euclidean distance between x and dD. It is well-known that any C > open 
set D with C 1,1 -characteristics (R,2,A) satisfies both the uniform interior ball condition and the 
uniform exterior ball condition: there exists ro = ro(i?2 ; A) < R2 such that for every x G D with 
o~8d(x) < ro and y £ ~R. d \D with <5ao(y) < ro, there are z^,^ G <9-D so that \x — z x \ = 5qd(x), 
\y - z y \ = 5 dD (y) and that B(x Q , r ) C D and 5(y , r ) CR d \D for z = ^ + r (x - z x )/\x - z x \ 
and y = z y + r (y - z y )/\y - z y \. 

In the remainder of this section, we fix such an ro and set To := $(ro/16). For any x G D with 
6d(x) < ro, let z x be a point on dD such that \z x — x\ = <5d(x) and n(z x ) := (x — z x )/\z x — x\. 

Lemma 5.1 Let kq G (0, 1) and a > 0. There exists a constant c = c(kq, R2, ro, a, (/>) > suc/i that 
for every (t,x) G (0,T ] x T> witt 5rj(x) < 3$ _1 (t) < r /4 and k G (0, 1), 

P x [Xg G B^, ko*" 1 ^))) > cy] * {5D ( x) \ (5.5) 
w/iere x := 2r x + \<&~ l (t)n(z x ) . 



Proof. Let < k\ < kq and assume first that 2 _4 Ki<J>~ 1 (i) < 5d(x) < 3$ _1 (t). As in the proof 
of Lemma 3.4, we get that, in this case, using the fact that \x — xq\ G [| fvQ*!? -1 ^), 6<I> _1 (i)], there 
exist constants Cj = ch(ki, ro, a) > 0, i = 1, 2, such that for all t < To, we have 



(x^eBh^jr 1 ^))) > cit($- 1 (t)) d J(x,x ) > c 2 >o. 



(5.6) 



By taking k% = ko, this shows that (5.5) holds for all a > in the case when 2 «o^ 1 (t) < 
S D {x) < 3$- x (i). 

So it suffices to consider the case that &d{x) < 2 -4 Ko < £ -1 (t)- We now show that there is some 
ao > 1 so that (5.5) holds for every a > ao and 5d(x) < 2 _4 ko < I >_1 (£)- For simplicity, we assume 
without loss of generality that xq = and let B := B(0, ko^ _1 (*))- Let U := D D B(z x , kq^ 1 ^)). 
By the strong Markov property of X D at the first exit time tjj from U and Lemma 3.2, there exists 
C3 = £3(0) > such that 

> F x (tu < at, X TU G ^(O^^ko^ 1 ^)) and \X? - X m \ < 2~ 1 k $~ 1 (*) for a e [m, t& + at]) 

> c 3 P x (t v < at and X r[7 G B(0, 2~ 1 K;o$ _1 (t))) . (5.7) 

Let x\ = z x + 4 Kon^a-)*!? (t) an d ^1 := -S( x i) 4~ 1 Ko$ _1 (^)). It follows from the boundary 
Harnack principle (Theorem 2.3) and (2.4) that there exist c& = Ck(R2, A, 7, 4>) > 0, k = 4, 5, such 
that for all t G (0,T ], 



\ (X^ G B(0, 2- 1 K $- 1 (t))) > c 4 P K1 {X TU G 5(0, 2- 1 Ko ^" 1 ^))) 

> c 5 F xl (x TBi G ^(O^- 1 ^^ 1 ^))' 
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By (2.10b, Lemma O and (T272]>-d2~4|>, we have 



X1 (X TBi G B(0,K $-\t)/2)) = E X 



JB(0,K 'S>- 1 (t)/2) 



J(X s ,y)dy 



> c 6 j(K <S>-\t)) 15(0, no*- 1 ®/!)} E x [t Bi ] 



> 



C7 



(At $- 1 (i)) d ^(^o^- 1 W) 



(«; $- 1 (i)) <i *(«o*- 1 (*)) = C7- 



Thus 



It follows from (5.3) that there exists eg > such that 



,(71/ > at) < (at)- L E x [Tu] < a- l cg 



*(8 D (x)) 



(5.8) 



Define ao = 2cg/(cg). We have by ( 5.7 >— (|5T8|) and the display above that for a > ao , 



G B) > c 3 (P,^ G 5(0, 2- 1 K $- 1 (t))) - P x (ra > at)) 



> c 3 (c 9 /2) 



$(M*)) 



(5.9) 



(5.6) and (5.9) show that (5.5) holds for every a > ao and for every x G D with Sjy(x) < (i). 
Now we deal with the case < a < ao and 5d(x) < 2~ 4 Ko$>~ 1 (t). If 8d{x) < 3&~ 1 (at/ao), 



we have from (5.5) for the case of a = ao that there exist cio = cio(fto> A, a, </>) > and 



en = ch(ko, R2, A, a, (p) > such that 

Pa^GB^MT 1 ^))) > V x (x° {at/ao) eB(xo, K o<S>- l (at/a )) 



$(5 D (x)) H8 D (x)) 
- c io\ — —, = Cn 



If 3$ _1 (ai/a ) < S D (x) < 2' i Ko^ 1 {t) (in this case 1 > k > 3 • 2 4 $~ 1 (a/a )), we get from 
(5.6) by taking k\ = $ (a/cto). The proof of the lemma is now complete. □ 



Proof of Theorem |1.5| (i). By Theorem |1.3[ i), it suffices to show that 

Px (t d >t)>c(lA !^MA 1 f or (t t x ) € (0, T] x D. (5.10) 



Assume (t,x) G (0, T] x Z). Since D satisfies the uniform interior ball condition with radius ro 
and < (T /T)t < T , we can choose a point as follows: if 5 D (x) < 3$- 1 ((T /r)t), let 
H = z x + {9/2)<S>- 1 ((T /T)t)n(z x ) so that 

5(4(3/2)<ir 1 ((To/T)i))c5(2 a; + ^^ 

and 5 D (z) > S^-^CTo/T)*) for every * G (3/2)cD- 1 ((T /r)t))). If «5 D (s) > SS-^CTo/r)*), 

choose ^ G B(x,5 D (x)) so that |x - f*| = (3/2)$" 1 ((T /T)^. Note that in this case, 

(3/2)d>- 1 ((T /T)t)) C 5(x, ^(x)) \ {x} 
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and 5 D (z) > $- l ((T /T)t) for every z G 2~ 1 <i>- 1 ((T /T)t)). We also define £* the same way. 

If <5 D (x) < 3$ -1 ((T /T)i), by Lemma [Tl] (with M = T/T when T > T ) and Lemma [57T| (with 



; (r D > t) > Cl F x (t d > (T /T)t) > Cl F x U D To/T)t G B(£, 2- 1 <S>- 1 ((T Q /T)t))) > c 2 



Hs D (x)) 



If 5 D {x) > 3$ _1 ((ro/r)t), by Lemma 4.1, Proposition 3.6 and ((O), 



x (td > t) > CiP, f Xg, /T)t G 2- i <D- i ((T /T)t)) 



ci 



B(€«,2-i#-i((2b/T)t)) 



p D ((T /T)t,x,u)du > c 3 > c 4 ( 1 A 



*(<5d(x)) 



1/2 



□ 



6 Large time heat kernel estimates 

In this section, we first give the proofs of Theorems |1.3[iii) and |1.5[iii) 



Proof of Theorem 1.3(iii). Since D is bounded, in view of (2.5), the transition semigroup 
{P t D , t > 0} of X D consists of Hilbert-Schmidt operators, and hence compact operators, in L 2 (D; dx). 
So Pf has discrete spectrum {e~ Xkt ;k > 1}, arranged in decreasing order and repeated accord- 
ing to their multiplicity. Let {4>k,k > 1} be the corresponding eigenfunctions with unit L 2 -norm 
H||l 2 (D) = 1) which forms an orthonormal basis for L 2 (D; dx). 
Clearly, for every k > 1 



D 



> l)4>k(x)dx < |Z?| 1/2 ||0 fc || L 2 (D) = \D\ 1 ' 2 . (6.1) 
By using the eigenfunction expansion of p£> we get 

f F x (t d > l)p D (t, x, y)F y {T D > 1) dxdy = V ( [ F x {t d > l)<f> k (x)dx\ . (6.2) 

JDXD \JD / 

Noting that is increasing and H/H^m) = YlkLiilo f( z ) ( t > k{z)dz) 2 , we have for all t > 0, 

/ F x (t d > l) PD (t,x,y)F y (T D > I) dxdy < e~ tXl [ F x (t d > l) 2 dx 
Jdxd Jd 



<e~ tM \D\. (6.3) 

On the other hand, by Theorem |1.3[ ii), Remark |l.l| ^iii) and (6.1) we have that there is a 
constant c\ > so that for every x G D, 



4> 1 (x) = e Xl / p D (l,x,y)4> 1 (y)dy 
Jd 

< Ci F x (t d > 1) / F y (r D > l)Mv)dy < ciIDI^PxCtd > 1). (6.4) 
Jd 
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It now follows from (6.2) that for every that for every t > 0, 



f F x (t d > l)p D (t,x,y)P y (T D > 1) dxdy > e" Ul f I F x {t d 
Jdxd \Jd 



> \)4>i{x)dx 



>e- tAl U c^Dr^Mxfdx 
For t > 3 and x,y £ D, we have that 



p D (t,x,y) 



p D (l,x,z)p D (t - 2, z,w)p D (l,w,y)dzdw. 



(6.5) 



(6.6) 



DxD 



By Theorem 1 1.3[ ii), Remark |l.l| ^iii), (2.5) and (6.3) we have that there are constants q > 0, i = 2, 3, 
so that for every t > 3 and x,y £ D, 

PD{t,x,y) < c 2 F x (t d > 1)P„(td > 1) / P 2 (r D > l)p (i - 2, z, w)P to (r D > l)cfedw 

JBxO 

< c 3 F x {r D > l)F y (r D > l)e~* Al . (6.7) 



By Theorem |L3^i), Theorem 3.7, the boundedness of I? and (6.5) we have that there are constants 
q > 0, i = 4, 5, so that for every i > 3 and x,y £ D, 

p D (t,x,y) > c 4 F x (t d > l)F y (T D > 1) f F z (r D > l) PD (t - 2,z,w)F w (t d > \)dzdw 

J DxD 

>c 5 F x (t d > l)F y (T D > l)e- tXl . 



This combined with (16.71) establishes Theorem 1.3 hi). 



□ 



Proof of Theorem |1.5| (iii). By Theorem |1.5[ i), it suffices to prove the theorem for T > 3. By 
( pUj ), JglQl ) and the boundedness of £>, P x .(r D > 1) X ^(^(x)) 1 / 2 . This and Theorem [l~3^iii) 
imply Theorem |1.5[iii). □ 



7 Green function estimates 



In this section, we use Theorem 1.5 to get sharp two-sided estimates on the Green functions of 
subordinate Brownian motions in bounded C 1,1 open sets. We first establish the following two 
lemmas. 

Lemma 7.1 For every r G (0, 1] and every open subset U ofM., 

2\ *(\x-y\) )-{ *(\x-y\y/*)\ *(\x-y\y/*) 

<lA r 2 $(^(x))V 2 $(^(y)) 1 / 2 
*(|x-y|) 
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Proof. The second inequality holds trivially. Without loss of generality, we assume 5u(x) < 5\j(y). 
If both r ^^ u ^ x l\ to and r ^ u ^ v )l L are less than 1 or if both are larger than one, 

/ r$(^(x))V 2 \ / r^ujy)) 1 / 2 ] r^ujx^/^ujy)) 1 / 2 

\ *<\x-v\yi*)\ *(\*-v\) 1/2 J *(\x-v\) 

So we only need to consider the case when ^ ^^ l/ jr < 1 < 7 ~^~~^ijy ■ Note that ${5u(y)) < 



$((5t/(x) + |x — 2/|). If <5f/(x) > |x -y\, then by (1.8), $(fo(y)) < $(26u(x)) < 4$(«Ji/(x)) and so 



" *(|ar-y|) " I <&(|* - y\)W ) ' 



When o"[/(x) < |x — y|, by (1.8) again, <&(5u{y)) < $(2\x - y|) < 4$(|x - y\) and so 



r^foKs)) 1 / 2 *^)) 1 / 2 < i a 2r 2 d>(^(x))V 2 / rd>(^(x))V 2 \ 
<*>(!* -y|) " *(|x-y|)V2 " ^ Sflx-yDV^' 

where the assumption r < 1 is used in the last inequality. This establishes the first inequality of 



(|7TJ). □ 

By condition (A), we have that for every T > 0, there exist Cr > 1 such that 



_ x f r x 1/(250 $-!(r) / rx l/(25 2 ) 



°t ( o) < / J, < Or ( „ ) forO<r<i?<T. (7.2) 



Moreover, for every M > 0, we have 

r$'(r)x$(r) forre(0,M] (7.3) 
(see the paragraph after |19| Lemma 1.3]). 
Lemma 7.2 Suppose T > and set 

h T (a,r) = a + $(r) [ (l A ^) - 1 rfu + (l A . (7.4) 

Then 



a I a I $(s) J 



for < r < <I>" 1 (r/2) and < a < (2" 1 A (2C T )" 252 )T ; w/iere C T is the constant in §L2§ and 
x+ ■= x V 0. 



Proof. For (a,r) with < a < $(r) < T/2, 

/*(r)/T w* -1 ^ -1 *^)) ' ' r < / 0(r)/T *- 1 («- 1 *(r))'* ' r 



h T (a,r) x a + a / _ + - = a+ - / \\ u du ' 
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By (ff^j), since $(r) < T/2, we have 



< c 2 = c x / u 2<5 i cZu < 
'1/2 



$-i($( r )) _ 1 



/*(r)/T 

Thus, for < a < <£(r) < T/2, we have 



it du < c\ I u 2S 2 du = C3 < oo. 



a . . / a\ a 

c 2 - < h T {a,r) < c 3 a + - < c 4 - 
r \ r / 



(7.5) 

r 

On the other hand, for (a, r) with $(r) < a < (2" 1 A {2Ct)~ 2&2 )T, using the change of variable 
u = <&(r)/<I>(s) and then applying integration by parts for the first integral below, we have 

t>$(r)/a 



h T (a,r) x a + $(r) / 

J^>(r)/a 



du 



+ a 



du 



a + 



ds + a 



a $(r" 



* ( T ) $'(s) , 3>(r) 

*-i(a) r 



+ a 



— + — 
*-i(o) s$ ( s ) r 



a + 



*-!(a) 



+ 



1(a) 



-ds + a 



(is. 



0-!(a) 



Since a < (2" 1 A (2Ct) _2<52 ) T, by ([T^ and the fact that is increasing, 

1 1 1 

> C4 



*-!(a) $-!(T) *-!(a) 



(7.6) 



(7.7) 



for some C4 > 0. Using (7.3) and (7.7) in the second integral in (7.6), we get that for (a,r) with 
*(r) < a < (2- 1 A (2C T )' 2&2 )T, 



h T (a,r) x a + ^ + / 



$(s) 



a + 



*- 1 (a) 



ds + a 



ds + a 



*-!(a) s 

1 1 



+ 



* _1 (o) 



$(s) 



(is. 



(7.8) 



Since $(s) is an increasing function, when < $(r) < a, we have 



*-!(a) 

while when <3>(r) > a > 0, 



+ 



$- 1 (a) 



$(s) , a 
— ^-ds < ,. 77 : + a 
s 2 $ _1 (a) 



+ 



*-!(a) 



$(s) 



1 , a 

^T«S = -, 



a a 
> -. 



a r 



This combined with (7.5) and (7.8) establishes the lemma. 



□ 



Recall that the Green function Gn(x,y) of X in D is defined as Gn(x,y) = PD(t,x,y)dt. 
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Theorem 7.3 Suppose that X is a subordinate Brownian motion with Levy exponent = 
<Ml£| 2 ) with 4> being a complete Bernstein function satisfying condition (A). Let D be a bounded 
C 1 ' 1 open subset ofM d with characteristics (R,2,A) anda(x,y) = $(5d(^)) 1 ^ 2 $(5_d(2/)) 1//2 , x,y € D. 

(i) There exists c\ > depending only on diam(D), R2, A,d and 4> such that for all d > 1 and 

(x,y) GDxD, 

(ii) There exists C2 > depending only on diam(D),B,2,A,d and <p such that for all d > 1 and 
(x,y) £ D x D, 

(iii) Let d = 1. T/ien /or (x, y) e D x D, 

GD{X ' V) ~\^\ A [*-Ha(x,v)) + ^7 J' 

(iv) Let d > 2. T/ien /or (x, y) € D x D, 

$(\x-y\) f lA a(x,y) 



x-y\ d \ §{\x-y\) 



Proof. Put T = (2 V (2C T ) 2,52 )$(diam(Z))), where C T is the constant in fl7!2| ). It follows from 
Theorem 1 1.5[iii) that 



/■OO 

/ PD(t,x,y)dt>c a(x,y). 
Jt 



(7.9) 



Using the boundedness of D, Remark |l.l[ iii) and (2.2), the results of Theorem |1.5[ i)-(ii) can 
be rewritten as follows: there exists c\ > such that for (t, x, y) G (0, T] x D x D, 



rl ( 1A ^™y /2 ( 1A ^(^M) i/2 ( ($ -i 



(t))- d A 



< PD(t,x,y) 



< ci 1 A 



$(Mx))\ 1/2 



1 A 



1/2 



($- 1 (t))- d A 



|x-y| d $(|x-y|) 

f ^ 

\x - y\ d ${\x - y\) J 



(7.10) 
(7.11) 
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By the change of variable u = y ^ and the fact that t — > <3? is increasing, we have 



' ^ ^(5 D (x)) \V 2 A A *(&(»)) 



o 



|x — y\ d 



mx-y\) 
\x — y\ d 

{ &(\x-y\) 

\x — y\ d 
--: I + 11, 



+ 



*(|*-y|)/T Ji 



u 



t 

-2 



1/2 



$-1(tQ 



A it 



-1 



x-y| d $(|x-y|) 

u$(M*)) 1/2 



df 



1 A 



{ H\x-y\y/2 J 



i 

*(|*-»|)/T 
oo 

-3 



-2 



I » - y| 



u 1 A 



$-i(7t-l$(|x-y|)) 
u^{8 D {x)Y^ 



1 A 



du 

ua(x,y) 



$(|a?-y|)V2 



1 A 



$(|x-y|) 
t*$(M2/)) 1/2 > \ 

^(1^-7/1)1/2 I 



f/it 



where in the fourth line of the display above, we used Lemma 7.1 
(i) The estimate on II is easy. 



l $(|x-y|) ( $(fo(x))V 2 \ / $(fo(y))V 2 \ 

" ' ^(k-2/l) 1/2 / V ^(k-y|) 1/2 y 

I $(|X- 7y|)V2 I I $(^-7/1)1/2 I 



2 |x-y| d 
$(|x-y|) 



<JJ < 



< 



X — 7/| d 

|x — y| rf 

|x — 7/| d 



du 



-1/2 



A 



<£(<M*)) 1/2 



u _i/ * A 



*(|x-y|)i/a 

2 / %(xT\ / $(M?/)) 1/2 \ 

^ $(|x- 7/|)i/2^ ^(la? - ^ 



g(|z-y|) L *(M*)) 1/2 > \ A *(My)) 1/a Y 



I a; - y\ 



$(|x-y|)i/2 ) 



$(|x-y|)i/2 

Now part (i) of the theorem follows from ( |7.10 ) and the lower bound of II in (7.13). 
(ii) We let 

a(x,y) 



u :-- 



$(|x-y|)' 



Clearly l/u > $(|x - y|)/$(diam(D)) > 2$(|x - y|)/T. By ( |7i2| |, 

a(x,y) f 1 k~y| d 



I < 



\x ~ y\ d U(\x-y\)/T ^-Hu-^dx - y\)Y 
a(x,y) 



7t 1 C?7i 



< c 3 

< C4 



x — y| d 
a(x,y) 



<5>- l m\x-y\)) 
<„, ,_„ , -,. V $-i(7i-i$(|x-y|)) 



Tt du 



|x - y| f 
x — y| d 



_2 i 

it 25 2 du 



(7.12) 



1/2 A ^Div)) 112 \ 



§{\x - y \yi 2 J 

du 



du 



(7.13) 



(7.14) 



(7.15) 
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Combining (7.9), ( |7.11 ), (7.13) and (7.15) we immediately get part (ii) of the theorem. 

(hi) Let h T (a,r) be defined as in (TO}. Since a < $(diam(L>)) < A {2C T )- 2S2 ))T, we have 



by (7.9)-(7.13) and Lemma 7.1 that Gr>{x,y) x hx{a{x, y), \x — y\). The assertion then follows 



from Lemma 17.2 



(iv) Note that since d > 2, we have by (7.2) that 

§(\x-y\) f 1 ___ 2 ( \x~ //I 



\x — y\ 



u 



*(|x-y|)/T 



< $(\x-y\) a(x,y) 



< 



C6- 



\X - J/|" 

H\x-y\) 
\x — y\ d 



1 A 



®{\x-y\)J Jo 

a(x,y) 
®(\x-y\ 



^(u-^dx-yl)) 



u d/(25 2 )-2 du 



ua(x,y) 
<f>(\x-y\) 



(7.16) 



Part (iv) of the theorem now follows from assertion (i) of the theorem, Lemma 7.1, (7.9), (7.11) 



(7.13) and (7.16). 



□ 



Corollary 7.4 Suppose that X is a one- dimensional subordinate Brownian motion with Levy expo- 
nent ^(C) = <KI£| 2 ) with 4> being a complete Bernstein function satisfying condition (A). Let D be 
a bounded C 1 ' 1 open subset o/M with characteristics (i?2,A) anda(x,y) = <^{5D{x)) l ^ 2 ^{5D{y)) 1 ^ 2 , 
x,y G D. 



(i) Suppose that for each T > there is a constant c\ = c\ (T, 0) > suc/i £/iai 



(is < ci 



<3?(r) 



/or r G (0,T]. 
a(x,y) 



F-y| V $(\ x -y\). 

(ii) Suppose that for every T > 0, there is a constant C2 = C2(T, <f>) > so that 

<f>(s) $(r) 

— ds < C2 for every r 6 (0,Tl. 

s^ r 

T/ien /or all (x, y) £ D x D, 

n ( n _ a(s,y) 
G °<™> * «-■(..(*,,,)) 

Remark 7.5 Recall that <5i,o~2 G (0, 1) are the constants in condition (A). 



(7.17) 



(7.18) 



(i) Condition ( 7.17[ ) is satisfied when 62 < 1/2. This is because for t G (0, T], we have by (1.7) 
and condition (A), 



$(r) 



r 



<D(s) 



(is 



r -2\ j /-T / s 2\ 5 2 1 

'-—ds<a 2 ^ ^(is = a 2 r- 252 



fs-2) S 2 



r z 



3 252 " 2 (is < 
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(ii) Condition (7.18) is satisfied when 5± > 1/2. This is because for t £ (0, T], we have by (1.7) 

(f>(r- 2 ) 1 



and condition (A), 



1 



<&(*) 



ds 



is' 2 ) s 2 



ds < a\ 



s 2 \ Sl 1 



ds = a ir - 2S2 I s 2Sl - 2 ds < -. 



Remark 7.6 Let cp(r) = r 1 / 2 /4>(r). Note that 

-l A i/2-& > ^{Xr)/ip(r) > a^X 1 ' 2 - 52 for A > 1 and r > R . 

Let 



<p*(X) := limsup <p(\r)/ip(r) and t/?*(A) := liminf (p(Xr) / (p(r) . 



Then 



oo > a^X 1/2 - Sl > <p*(\) > <^(A) > a 2 l X 1,2 - &2 > for A > 1. 
The upper and lower Matuszewska indices can be computed as 

a(ip) := lim (log <p.( A)) /(log A), /%) := lim (log cp* (A)) /(log A). 

A— >oo A— >oo 

Thus we have 

oo > a^A 1/2 - 51 > y>*(A) > A^ > X a ^ > ^(A) > a 2 1 X 1 / 2 - &2 > for A > 1 

(see [HI page 69-71]). Let <p{r) = j£_i /2 <p(t)/tdt, r > T~ 2 . With the change of variable s = r 1 / 2 , 
we see that ( 7.17| ) is equivalent to 



ip(r) 



r 
oo 



0(t)tV2 



dt = 2 



*00 



-1/2 S z 



< ci $(r^ 1 / 2 )r 1 / 2 = <p(r) for r > T~ 2 . (7.19) 



Let £>(r) = J ' ip(t)/tdt. (7.18) is equivalent to 



<p(r) 



-1/2 ( 

2 / -^ds < c 2 $(r- 1 / 2 )r 1 / 2 = <p(r) for every r > T~ 2 . (7.20) 
in s 



0(t)tV2 



dt 



In fact, by j3j Corollaries 2.6.2 and 2.6.4], ip(r) x ty?(r) f° r every r > T- 2 if and only if /%) > 0, 
and </?(r) x ¥>(r) for every r > T~ 2 if and only if a(ip) < 0. 

By following the same proof in |21|, Section 6], one can construct <f> whose upper and lower 
Matuszewska indices are a{4>) = 3/4 and /?(</>) = 1/4 so that a(<p) = 1/4 and /3(<p) = —1/4. For 
such 4>, neither (7.17) nor (7.18|) hold. 



Proof of Corollary |7^4| Put T = (2 V (2Cr) 2<52 )$(diam(D)), where diam(£>) is the diameter of 
D and Ct is the constant in (7.2). For notational simplicity, we let a = a(x,y) and r = \x — y\. 
Recall that uq 



a/$(r) is defined by (7.14). In view of Theorem |7.3[ iii), we only need to consider 
^0 > 1) which we will assume from now on in this proof. 
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(i) Since a < $(diam(L>)) < (2C T )~ 252 T, by (|7Jj) and (ill]) we have 



+ 



-ds x a 



1 



i \ r*' 1 ^ $( s ) , 



= a/ -5-+/ A^dsK Ards<cx^- 



Combining this with Theorem |7.3[i), (iii) and Lemma 7.1 establishes part (i) of the corollary. 



(ii) Since a < *(diam(D)) < T/2, by (7.18) we have that for u > 1 



a 

*-!(a)' 



Thus by Theorem |7.3[iii) , 



G D (x,y) 



$-!(a) 



when uq > 1. 



Since <3?(r) is increasing in r, the above together with Theorem 7.3 iii) implies that 



o a Jo/r if a < $(r) 

r &-Ha) \a/Q- l (a) ifo>*(r). 



□ 



We next give an example of a one-dimensional subordinate Brownian motion with Levy exponent 
0(|£| 2 ) that satisfies condition (A) but its associated function <3?(r) = l/(p(r~ 2 ) satisfies neither 



condition (7.17) nor (7.18). We can get its explicit Green function estimates on bounded C > open 



sets by using Theorem |7.3[iii) but not by Corollary 7.4 



Example 7.7 By the discussion at the end of the Introduction, we know that for any p E R, the 
function 



J {r + ty 2 t l ' 2 {\n{t + l)fdt\ , r>0 



is a complete Bernstein function. Moreover, 0(r) x r 1 / 2 (logr) p when r > 2. When p E [—1/2, 0) U 
(0, 1/2], (j)(r) = r 1 / 2 (log(l + r)) p is an explicit example of a complete Bernstein function such that 
4>(r) x r 1 / 2 (logr) p when r > 2. Let X be a one-dimensional subordinate Brownian motion with 
Levy exponent </>(|£| 2 )- Note that 



$(r) = l/</>(l/r 2 ) x r(log(l/r))~ p for < r < 1/2. 



(7.21) 



Let W(x) be the Lambert-W function, that is, W(x) is the unique solution to x = W(x)e w( - X \ 
Suppose s = r(log(l/r))~ p . Then 



y --p-^s-VP = r- l l' p \og({l/r) l / p ) =: 



xe 



So logr- 1 ^ 



= W(y) = W(p 1 s 1 ' p ), which implies that r = exp(— pW(p 1 s 1 ^ p ))- So in view 
of (1.8), there is a constant cq E (0, 1) so that 

Co exp(-pJ^(?T V 1 ^)) < $ _1 (s) < eg 1 exp(-pVF(jrV 1/p )) for s E (0,$(l/2)]. (7.22) 
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Recall that 



(log(l/ S ))- 



-ds 



1 



and 



p-1 
(logil/s))- 1 



(log^/s)) 1 p + c, when p / 1 



ds = -log(log(l/s)) + c. 



Suppose < r < $ -1 (a) < 1/2. Then if p ^ 1, 



ds 



<S>- 1 (a) 



(log(l/ S ))-P 



1 ((log(l/$- 1 (o))) 1 ^-(log(l/r)) 1 ^) 



(7.23) 



and, if p = 1, 



(is 



(log(l/ S ))- 



ds = loglog(l/r) — loglog(l/$ 1 (a)) 
log(logr/log$ _i (a)). (7.24) 



s 

,-1/ 



Let Z) be a bounded C ' open set in M. We further assume that 

$ -1 (diani(.D)) V diam(D) < c /2, 

where cq G (0, 1) is the constant in ( |7.22| ). 

Let a(x,y) = ^(5 D (x)) 1 / 2 ^(5 D (y)) 1 / 2 . We have by Theorem Qui), ff^Ij ) and ([733]) -([734) 
that for p / 1, 



|x-y| 1$ x (a(a;,y)) \J\ x - y \ 



\x - y\ 



A lo. 



<&- l {a{x,y)) 



(logil/Q-Hajx^)))) 1 p - (log(l/|x - y\)) 
p — 1 



+ 



while for p = 1, 

G D (x,y) x ^4 A | ! lo, 



1 



+ log+ (log|x-y|/log$ 1 (a(x,y))) 



where log + x := V logx. It is elementary to check that for < u, r < cq/2 and cqu < v < Cq 1 n, 
1V P / (log(l/«)) 1 ^-(log(l/r)) 1 - p V 



log- 



it 



p — 1 



when p / 1 and 



log—] + log + (logr/logu)) x ( log - ] + log + (logr/logv)) . 
uj \ vj 
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Thus we have from the four displays above together with (7.22) that for p 7^ 1, 

\x-y\ \ \ p-i 

while for p = 1, 

G D (x,y) x ^^A(w(a(x,y))- 1 + log+(iy(a(x,y))- 1 log(l/|x-y|))). 



Acknowledgement. The main results of this paper, in particular Theorems |1 . 3| and Corollary [L6J 
were reported at the The Sixth International Conference on Stochastic Analysis and Its Applications 
held at Bedlewo, Poland, from September 9 to 14, 2012. At the same meeting, K. Bogdan, T. 
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in terms of surviving probabilities and global heat kernel p(t, c\x— y\) as in our Theorem |1.3[ i)(ii) for 
the Dirichlet heat kernels of a similar class of purely discontinuous subordinate Brownian motions 
considered in this paper but only for bounded C 1,1 open sets. 
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